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Abstract 

In order to investigate the charge symmetry breaking (CSB) in the short range part of the 
nuclear force, we calculate the difference of the masses of the neutron and the proton, AM, the 
difference of the scattering lengths of the p-p and n-n scatterings, Aa, and the difference of the 
analyzing power of the proton and the neutron in the n-p scattering, AA{9), by a quark model. 
In the present model the sources of CSB are the mass difference of the up and down quarks and 
the electromagnetic interaction. We investigate how much each of them contributes to AM, Aa 
and A^(^). It is found that the contribution of CSB of the short range part in the nuclear force 
is large enough to explain the observed AA{6), while Aa is rather underestimated. 

PACS numbers: 12.39.Jh, 13.75.Cs, 13.88. +e, 13.40.Ks 
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I. INTRODUCTION 



The charge symmetry is the invariance under the charge-reflection, i.e., the reflection 
about the 1-2 plane in the isospin space. If this were an exact symmetry, the masses of the 
proton and the neutron would be the same, as well as the binding energies of the mirror 
nuclei or the scattering lengths of the p-p and n-n scatterings. The charge symmetry holds 
only approximately in the real world. There are small but non-zero differences such as 

AM = M„ - Mp = 1.29 [MeV] and Aa = app - ann = 1-5 [fm] . (1) 

These differences are manifestation of the charge symmetry breaking (CSB). 

CSB appears also in spin- dependent observables. For example, the p-n system is the 



mirror of ii-p, where p(n) is a polarized nucleon. There was 
analyzing powers of p and n in the medium energy scattering 



bund small difference in the 



AA{9) = Ani9) - Apie) . (2) 

The study of AA{6) is important because there is no Coulomb interaction between n and p. 

It is important to understand CSB from the quantum chromodynamics (QCD) viewpoint 
j^. From QCD we find that CSB has two origins: (i) the difference of the masses of the up 
and down quarks and (ii) the electromagnetic interaction. Thus the study of CSB phenomena 
can be a good probe to examine the behavior of the quarks and gluons in the low-energy 
region. The ultimate goal of the CSB study may be understanding their effects on hadron 
spectra and hadronic interactions directly from QCD, by, e.g., lattice QCD simulation. As 
the direct approach is not available up to now, however, indirect approaches have been taken 
for the CSB study. 

An often used approach to CSB is based on the meson exchange picture of the nuclear 
force. It was suggested that CSB of the nuclear force is generated by mixings of / = and 
1 = 1 mesons such as p-u mixing |^. A model based on such a picture was reported to 
explain Aa well. P| But it was also pointed out that the effect of the p-u mixing to CSB 
may be suppressed by the off-shell effect of the p-uj mixing. [6] Thus, this problem is still 
open Q. A class IV interaction |8|] is also generated by the neutron-proton mass difference 
in the one-pion-exchange interaction. It was pointed out that the effects of OPE and p-u 
mixing explain AA{6) fairly well. 
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On the other hand, CSB appearing in the short-range part should be investigated by 
introducing subnucleonic degrees of freedom. One of the pioneering works to apply a quark 
model to CSB is found in Ref. where the isovector mass shifts of isospin multiplets and 
the isospin-mixing matrix elements in IsOd-shell nuclei are investigated by using the quark 
cluster model (QCM) U , 13, 13 13 • K was concluded that the u-d quark constituent 
mass difference produces significant effects, which may explain the observed Okamoto-Nolen- 
Schiffer anomaly well. 

In the present work, we investigate CSB in AM, Aa and AA{6) by employing essentially 
the same model for all these three observables: a quark potential model for AM and QCM 
for Aa and AA{6). The CSB sources are taken to be (a) the difference of the masses of 
the up and down constituent quarks and (b) the electromagnetic interaction between the 
constituent quarks. Our aim is to estimate the effect of CSB sources (a) and (b) on nuclear 
force by investigating the above three observables simultaneously. 

Chemtob and Yang 3l (CY) calculated Aa using QCM, s ugg esting that the quark mass 
difference contributes to Aa significantly. Later, Brauer et al.jlal3 studied Aa and AA{9) 
using QCM and concluded that the effects of CSB sources (a) and (b) are too small to 
explain the observed value. However, their calculation of AA{9) suffers from a wrongly 
chosen factor, from omitting the symmetric spin-orbit term and from inconsistent use of the 
operators and wave functions (See sec II V|) . 

In the present paper, we extend CY's and Brauer's works in order to obtain more inte- 
grated knowledge on CSB. We investigate CSB in AM, Aa and AA(£) simultaneously. Also, 



we introduce the Instanton Induced Interaction (III) |2ij|, |2l|, |22 , which comes 

from the nonperturbative effects of QCD and explains the t] — rj' mass splitting . Since III 
does not break the charge symmetry, its role in this study is mainly to make the effective 
strength of the one-gluon exchange interaction smaller. The strength becomes reasonably 
small, which fits to the picture that this term represents the perturbative effect of the gluons 
(See sec II V|l . Moreover, we include the symmetric spin-orbit term in the analysis of AA{6), 
whose effect is as large as the antisymmetric one. Furthermore, we solve QCM to obtain the 
relative wave function and use it to evaluate the matrix elements of Aa and AA{9). 

In section |nj we show the Hamiltonian for quarks and the CSB sources. In section ITTTl 
we explain the detail of the calculations of AM, Aa and AA{9). Results are discussed in 
section HVl Summary is given in section IVl 
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II. HAMILTONIAN 



We employ the constituent quark model with quark masses of order m ~ 300 [Me V] in 
this study. The Hamiltonian is given by 

H = K + V (3) 

K is the quark kinetic energy and considered as semirelativistic in calculation of AM (See 
sec 1111 A|) and as non-relativistic in calculations of Aa and AA{6) (See sec 1111 B|) in this 
study. The quark-quark interactions are represented by a static potential, which consists of 
the confinement (CF), the one-gluon-exchange (OGE) 2^, the electromagnetic (EM) and 
the instanton induced (111) interactions. 

V = Vconf + VoGE + Vem + Vm (4) 

i<j 

~ [tTtI 2 "I 2 ~^ ' 7^ ~ 

2rfj mi rrij minij 2 



4r?- mf 2 
Vem = Ee.e,«,„{ ^ " ^2^ + 2^ + 3^/^ " ^^^^^-^ 



mi mj mimj 



Arfj mf mj 2 



Vm = Vi'^ + ^A, . A, + ^A. ■ X,a, ■ a,)6{n,) 

i<j 

Aj is the color SU(3) Gell-Mann matrix and Cj is the quark electric charge in units of the 
proton charge e. In this study it is assumed that the confinement potential does not break the 
charge symmetry. This is a natural assumption based on the confining potential obtained, 
for instance, from lattice QCD calculation. Yet there may exist velocity dependent terms 
associated with confinement which break the charge symmetry. We do not consider such 
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terms in this study. Taking the Breit-Fermi interaction naively, non-Gahlei invariant terms 

appear in the LS terms. But we consider only the Galilei invariant terms such as LS term in 

Eqs. (jniEj). It should be noted that the Instanton Induced Interaction (III) is effective only 

on the flavor singlet (iso-singlet) quark-quark state. In other words, it works only on a pair 

of up and down quarks. Thus III does not break the charge symmetry. 

In this Hamiltonian the terms including the quark mass and the electric charge may break 

the charge symmetry. In order to show the CSB terms explicitly we rewrite the quark mass 

and the electric charge in terms of the isospin operator. 

md + rriu rud - rUu m 

rrii = Tn 

2 2 ^ 



m(l 



Am 



2m ^ 

-ii) 



m{l-ern (9) 



1 

6 



(10) 



where 



md + mu . 
m = Am = md — mu 

Am 



(11) 



2m 

Using the typical constituent quark mass rh ~ 300 MeV and the up and down quark mass 
difference Am ^ 6 MeV, e ~ ^xsoo ~ small as the electromagnetic coupling 

constant, ae.m. — 1/137. So we divide the Hamiltonian into the charge symmetric part H 
and the charge symmetry breaking part AHcsb, and treat AHcsb perturbatively. 
The CSB part of the Hamiltonian is given to the leading order in e and ae.m. by 

AVcsB = AYgs^B^^ + AV™ (12) 
AVgs^l = E(A.-A,)%{ -|j(r«+ri^'^)(l + ^a.-c?,)5(rl,) 

i<j 



a.. 



Lir(c?.-^.)(rf -r?)} (13) 



Am'^rfj 



AV™ = E ^^^^«e.m.{ - -^(1 + h ■ -Mr^.) 

i<j -"-^ 'ij -J 

^Lij-(cr. + c?,)} (14) 



We ignore the second order terms C(e^, , eae.m.)- The CSB terms from the tensor 
interaction are excluded because the tensor interactions between quarks are smalL But 
we consider them when solving the charge symmetric equation for the unperturbated wave 
function. 

The Hamiltonian has 5 parameters, as,fh,a,VQ^^ and Am. The parameters are deter- 
mined so as to reproduce the single baryon properties and the results are shown in section HVl 



III. CALCULATIONS 

In this section we present the formulas of the neutron-proton mass difference, AM, the 
difference of the scattering lengths of the p-p and n-n scattering, Aa, and the difference of 
the analyzing power of the neutron and the proton of the n-p scattering, AA{6). 

A. The proton-neutron mass difference AM 

The differences of the mass of the isodoublet hadrons were evaluated in the constituent 
quark model by Isgur [2^. We also evaluate the neutron-proton mass difference in order 
to determine the mass difference of the up and down constituent quarks. Our approach is 
different in the following two points. First, we consider the semi-relativistic kinetic energy 
term, 
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K = E\/^'+P' (15) 

i 

Eq. (fT3j) can be divided into the charge symmetric part and the charge symmetry breaking 
part, 

K = K + AKcsB (16) 
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K = J2^m^+pj (17) 



_ 2 
m 



AKcsB = - E /'^ (18) 



Eq. ()15|) contains the kinetic energy of the center of mass coordinate, which must be sub- 
tracted in order to calculate the baryon mass. For the semirelativistic kinematics, the 



center-of-mass energy can not be treated exactly. Therefore we use the following approxi- 
mation, 



Mn = - PI) 



(H) - g (19) 



(Pel 

better than that in (^). Then the nucleon mass can be written in terms of H and AH,^^^ 



The relativistic effect is partially included as the convergence of the expansion in ^"-^ is 



as 



Mn = (H) - ^ + (AHcsb)(1 + ^) (20) 



where 



H = H + AHcsB (21) 
H = K + V (22) 
AHcsB = AKcsB + AVcsB (23) 

H is the charge symmetric part of the Hamiltonian and AHcsb contains Eqs. ()18p and 
fll2|) . The first two terms of Eq. (PUI) give the average mass of the nucleon and the third 
term contributes to AM. The up-down quark mass difference Am is determined so as to 
reproduce AM by using Eq. (pUj) . 

The second difference from the Isgur's work is that the Instanton Induced Interaction (III) 
is considered in this study. Ill has the contact spin-spin interaction and contributes to the 
difference of the masses of the Nucleon and A(1232) just like the color magnetic interaction. 
We choose the coupling constant of the OGE, a^, and the III, , so as to reproduce the 
nucleon-A mass difference in total. So as becomes smaller effectively by considering III. 



B. CSB in the N-N scattering 



In the calculation o 
cluster model (QCM) 



le scatterin g le ngths and analyzing powers, we employ the quark 
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atterm g le i 



which describes two-nucleon systems in terms of 
their quark coordinates. The scattering wave functions, which are used as the unperturbated 
states, are calculated by solving the resonating group method (RGM) equation. By mainly 
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technical reasons the kinetic energy term is treated purely in the non-relativistic way, i.e. 
the semirelativistic kinematics is not taken into account contrary to the case of single baryon 
mass. This approximation can be justified because the relativistic effect on the kinetic energy 
term is smaller for the motion of the two baryons. Then the kinetic energy is given as 

K-X^Ki-Kc (24) 

i 

K. = K + ^) (25) 

Kg = ^ (26) 
2Mc ^ ' 



where 



MG^Y.mi Po^YlPi (27) 

i i 

The RGM equation for the baryon A and baryon B is as follows, 

/ MMM^bKU - E)AlMU)MiB)x{RAEmAd^E = (28) 

(f)A{B) and Ca{b) is the internal wave function and coordinates of the baryon A(B). Rab is the 
relative coordinates of the baryon A and B. The parameter b is the gaussian size parameter, 
which represents a nucleon size. A is the antisymmetrization operator for six quarks and is 
written as follows. 

Yl Pij (30) 
ieAjeB 

In the end, the following equation is obtained 



2AtAB 2/iAB 

- EN^^''\R,R'))xiR')=0 (31) 

where Pab is the momentum operator of the relative motion of the baryons A and B, and 

k2 

= Ma + Mb + —— (32) 

J- ^ -L + _L (33) 

AtAB Ma Mb ^ ^ 
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Ma(b) = (34) 
jeA(B) 

111 

^ = ^ + ^ (35 

/UAB Ma Mb 

Ma(b) : observed mass of the baryon A(B) (36) 

It should be noted here that Ma(b) and Ma(b) may not agree with each other completely. 
We take rrii = 313 [MeV] in our calculation so that the difference is small, but for the charge 
symmetry breaking we assume that /^ab = /^ab- The observed masses of the proton and 
neutron are given by 

Ma = M(l - eNrf ^) (37) 
AM 

CN = ^ 38 
2M ^ ^ 

M = Mid:^ = 939 MeV (39) 
AM = Mn - Mp = 1.29 MeV (40) 

Therefore we may rewrite the kinetic energy terms as 



2yUAB 2/iAB 2yU 

and the energy in Eq. as 



P''^ ^\l + ltl±fL,^) (41) 



because 



The RGM kernels vif\ N^^^\ K^^^\ V^^^^ are defined by 

Vi?{R) = J deAdeBdRAB0A(a)0B(eB) 



(42) 



^(A) (B) ,2 ^(A) (B) 

E = 2M(1 - "-^^^e.) + |(1 + "-^^^e,) 

1.2 L.2 AA) , (B) 

= 2M + — + -2M + —)(- — ^^eN 

= E + AEcsB (43) 



1 ^(A) (B) 

^ V' ^n) (44) 



2/iAB 2/i 2 

A ^ I (45) 



V,,(5(R-RAB)0A(a)0B(eB) 



'^^ac?^bc?Rab</'a(^a)0b(^b)5(R' - Rab) 



(46) 



(l \ 

K 



^'[5(R-RAB)0A(a)0B(eB)] 

= i?) + AiTcsB {R', R) (47) 

^V(^^)(i?',i?) + AVbsB(i?',i?) y 

K and V are given by Eqs. and Q and can be divided into the charge symmetric part 
K,V and the charge symmetry breaking part AKcsbj '^Vcsb- Therefore RGM kernels are 
divided into the charge symmetric part K^^^\V^^^'^ and the charge symmetry breaking 
part AKg^\Avg^\ 

In order to treat the CSB part perturbatively, we employ the distorted wave Born ap- 
proximation (DWBA) in this study, we solve the following equation to obtain the distorted 



wave. 



1^ - 



dB!{K^-^^\R, R') + R') 
EiV(^^)(/?,i?'))XdistW = 



(48) 



The direct kernel V^^\R) comes from the electromagnetic interaction of quarks and corre- 
sponds to the electromagnetic interaction of baryons. We are interested in effects of CSB 
at the quark level, not at the hadron level. So we ignore the direct kernel. But we consider 
the exchange kernel of the electromagnetic interaction of quarks. Using the distorted wave 
Xdis(-R), we estimate the following CSB parts. 



(CSB part) 



^AB ^ 



-(A) 



(B) 



-eN)Xdist(-R) 



2/2 ' 2 
dIi[AKgl\R, R') + AVg^\R, R') 

AEcsBiV(^^)(i?,/?')]Xdist(i?') 



(49) 
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C. CSB in the analyzing power 



There is a special CSB interaction in the neutron-proton system, which is called the class 
IV interaction, according to the classification by Henley and Miller jsj. 

Fiv (X {T^-Ti){aA-ffB) (50) 
or 

(rA X tb)z{^a X ae) (51) 

one sees that the class IV interaction mixes spin-singlet states and spin-triplet states. The 
spin singlet-triplet mixing induces asymmetries of spin polarization observables such as 
the analyzing power. At the level of the quark-quark interaction, CSB in the spin-orbit 
interactions is given as [See Eqs. (OZI)] 

■\/LS _ -irOGE I -irOGE I -irEM /cr,\ 
V CSB — ^ qSLS + V qALS + ^ qSLS V^^/l 

Vq^Ll = - E(A. • A.)|^^[(^. + ^.)(rf + 4^')] (53) 



V?aTs = - E(A. ■ - ^.)(rr - 4'^)] (54) 

i<j ^^'"^ 'ij 



Kls = - E lli^T? [(-^ + + -i^')] (55) 

i<j 'ij 

The first two terms of Eq. come from the one-gluon-exchange interaction and the third 
term from the electromagnetic interaction of quarks. It should be noted that the symmetric 
spin-orbit interaction of quarks (qSLS) induces the class IV interaction of baryons as well 
as the antisymmetric one (qALS). Brauer et al. calculated AA{6) using a similar model 
without including the qSLS terms Q]- They concluded that the contribution of quarks to 
AA{6) is very small. But we will see that the contribution of quark spin-orbit interactions, 
Eq. (|^. to AA{9) is large enough to reproduce the observed AA{9). 
Using DWBA, we calculate the following matrix elements for J=L<3, 

ATcsB = (^LjlV^lBrij) (56) 

Then the total T-matrix is given as follows, 

T = Tcs + ATcsB (57) 
11 



Tcs is obtained by solving the RGM equation. T is the regarded as a matrix based on the 
spin states and the analyzing power is given by 

_ TrfTt^NT] 

- Tr[TtT] ^^^^ 
Then AA{6) is given in terms of Tcs and ATcsb 



AA{e) = A„(0) - Ap{e) 

^ 2ReTrfe(an-ap) ATcsb] 
TrfeTcs] 

We show the explicit forms of the T-matrix and of AA{6) in Appendix A. 



(59) 



IV. RESULTS 



The parameters in our calculation are determined so as to reproduce the single nucleon 
property. In order to show explicitly how much the contribution of the Instanton Induced 
Interaction (III) to the Nucleon- A splitting is, we introduce a new parameter Pm, which 
denotes the ratio of the contribution of III to the whole Nucleon-A splitting. For example, 
when Pni = 0.4 the contribution of III to the Nucleon-A splitting is 40% of the whole one. 
Vq is determined so as to reproduce the r] and i]' mass splitting. Our analysis shows that 
Fill ~ 0.4-0.5 gives the right rj-ri' splitting. Here we try two values Fm = 0.4 and 0.5. Using 
the nucleon mass formula Eq. (j^Uj) . we obtain Am for each Fni. The results are given in 
Table lU The parameter b is the gaussian size parameter for the internal wave function of 
the nucleon, which represents the nucleon size. 



TABLE I: Parameters 





Piii 


Am 


m [MeV] 


b [fm] 


as 


a [ MeV/fm] 


V^^^ [MeV fm^] 


A 


0.4 


7.3 


313 


0.6 


0.91 


44.29 


-177.2 


B 


0.5 


5.2 


313 


0.6 


0.76 


40.34 


-221.5 



Another possible source of the N-A splitting is contribution of pion cloud around the 
)aryon. For instance, the cloudy bag model predicts the N-A splitting of about 100 MeV 
2^. This effect may reduce the roles of OGE and III, but it is not taken into account in 
this approach. 
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By increasing Pni, we reduce as accordingly so that the N-A mass difference is fixed. 
For Pni = 0.4, as becomes 0.91, while as = 1.52 is necessary to reproduce the N-A mass 
difference only by OGE. In order to show the effect of the Instanton Induced Interaction to 
AM, we show contribution of each term to AM in Table ITTl for various Pm. The Kin, OGE 
and EM represent the contributions of the kinetic energy, the one-gluon exchange interaction 
and the electromagnetic interaction to AM. It should be noted that when Pm = we can 
not reproduce the AM because OGE gives large contribution, which goes to the opposite 
direction. This shows the essential role of the III, which reduces the OGE strength. 



TABLE II: Contributions to AM for A to =6 MeV 

Pm Kin OGE EM M^ - Mp 

4.72 -5.54 -0.41 -1.23 

0.1 4.72 -5.54 -0.41 -0.67 

0.2 4.72 -4.99 -0.41 -0.12 

0.3 4.72 -4.43 -0.41 0.44 

0.4 4.72 -3.88 -0.41 0.99 

0.5 4.72 -2.77 -0.41 1.54 



It is also found that the calculation of "Strong hyperfine" for "p-n" in Table I of Ref . |2Sl 
is different from our calculation even if we use the same potential. This is because Isgur 
considers distortion of the quark wave function from the u-d quark mass difference. However, 
to be consistent the distortion of the wave function should not contribute to the energy in 
the first order of the perturbation theory. ^ The contribution of the "Strong hyperfine" to 
AM should be ^5^^ instead of -^S^?^ in Ref. 1281, where 5 is the nucleon-A mass splitting. 

Next we calculate Aa using the parameters in Table HI The results are shown in Table 
mil a(f) and Aa(Ar) is the average and the difference of the scattering lengths (effective 
ranges) of the p-p and n-n scatterings. 



2 



2 

^ Chemtob and Yang also point out the mismatch with Isgur in their pape [l^ . 



Aa = app - ann (60) 

Ar = Tpp - Tnn (61) 
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Our results, Aa =0.79 and 0.52 [fm] for Pni =0.4 and 0.5, are somewhat smaller than the 
observed value ~ 1.5 [fm]. We, however, point out that Aa is sensitive to the parameters 
because it is given by a cancellation of positive and negative terms. 



TABLE III: Scattering length 





Piii 


Am [MeV] 


a [fm] 


Aa[fm] 


r[fm] 


Ar [fm] 


A 


0.4 


7.3 


-17.9 


0.79 


2.42 


-0.39 


B 


0.5 


5.2 


-17.9 


0.52 


2.46 


-0.25 


Exp [3] 






-18.1±0.5 


1.5±0.5 


2.80±0.12 


0.10±0.12 


B [18j 




5.0 


20.07 


0.46 






CY [17] 




6.0 




2~3.5 







In Table CYl we show each contribution of CSB terms Eq. ^ to Aa. NMD, Kin, OGE 
and EM are contributions of the first term of Eq. (jl^ . the quark kinetic energy (including 
the AEcsB term), the one-gluon exchange interaction and the electromagnetic interaction, 
respectively. 

TABLE IV: The contributions to Aa of CSB terms [fm] 





NMD 


Kin 


OGE 


EM 


A 


0.3 


-2.6 


2.9 


0.2 


B 


0.3 


-1.7 


1.7 


0.2 



We estimate Aa in our formulation for the parameters of Ref. tlS, (B) andQ (CY). The 
contributions of Kin and OGE should be given by 



AaKin OC —^T-r = A^Kin (62) 

a^Am 

AaoGE OC _ = A60GE (63) 

In Table V we show A^Kin and A6oge for the parameters of B and CY. Using the values of 
Table HVl and El we find 

4.6 5.2 

AaKin + AaoGElB = -2.6 x — + 2.9 x — 

o.U (.1 

= -1.5 + 2.0 = 0.5 (64) 
14 



AoKin + AaoGEic = -2.6 X — + 2.9 X — 

= -1.8 + 3.4 = 1.6 (65) 

These estimates suggest that our resuhs may become larger by the changing the parameters. 
As Afexin is larger than A6oge in our parameter choice, the cancellation of Aaxin and Aooge 
is stronger than the other cases. On the other hand Ar is too large and has the wrong 
sign. More investigation should be done for Ar, which reflects not only the strength of the 
interaction but also its radial dependence. 



TABLE V: A^Kin and A^oGE 




A6Kin [MeV] 


A&oGE [Mev] 


A 


8.0 


7.7 


B 


4.6 


5.2 


CY 


5.6 


8.9 



Finally we calculated AA{6) at two energy points, taking Pni = 0.4. The results at En = 
183 and 477 [MeV] are shown in Fig^ □ and H The results at =183 and 477 MeV are 
lar ge e nough to reproduce the data lly], which disagrees with the conclusion of Brauer et 
al. [19| . The difference mainly comes from two points. The first point is that they consider 
only the antisymmetric spin-orbit interaction of quarks (qALS) not the symmetric spin-orbit 
interaction of quarks (sSLS). The factor of qSLS is three times as large as that of qALS (See 
Eqs. fj53ll54j) ). The remaining discrepancy might be attributed to their erroneous choice of 
the unit of 71 in the formula Eq. (3.6) in their paper ^|. We convert their value of 71 in 
radian into that in degrees and obtain AA{6 = 96°) = 5.4 x 10~^, which is of the same order 
as our estimate. Our result at En = 477 [MeV] is too large. It is not surprising since we fit 
the phase shift of the N-N scattering up to En = 400 [MeV] and we may not apply QCM at 
higher energy and we need higher partial waves. 

Fig. El and H show the contributions of Pi\f\'-^ Pi) , (^Dalf j^Z^z) and (^FglfpFg) to 
AA{6). It is found that the contribution of {^Pi\T\^Pi) is dominant in the observed 6 
region. But the other mixings of partial wave become important for the other 6 region. 

We also investigate each contribution of the one-gluon exchange interaction and the elec- 
tromagnetic interaction. (Fig. El and Ej) It is found that the contribution of OGE depends on 
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FIG. 1: A^(^) at En = 183 [MeV]. 
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FIG. 2: AA{e) at En = 477 MeV. 

the incident energy much strongly than that of the electromagnetic interaction does. This is 
because the dominant contribution of the EM interaction is the direct interaction while OGE 
interaction contributes as the exchange interaction. Therefore their energy dependences are 
different from each other, which may be studied by future experiment at various energy 
points. 
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FIG. 3: The contribution of each partial wave mixing at En = 183 MeV. 
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FIG. 4: The contribution of each partial wave mixing at En = 477 MeV. 
V. CONCLUSION 

We have calculated the difference of the masses of the neutron and the proton, AM, the 
difference of the scattering lengths of the p-p and n-n scatterings, Aa, and the difference 
of the analyzing power of the proton and the neutron in the n-p scattering, AA{9), using 
the quark cluster model. In the calculation of AM, we treated the kinetic energy in the 
semirelativistic way and introduce the Instanton Induced Interaction (III). We have found 
that the contribution of the one-gluon-exchange interaction (OGE) is suppressed by the 
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FIG. 5: The contribution of OGE and EM at En = 183 MeV. 
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FIG. 6: The contribution of OGE and EM at En = 477 MeV. 

introduction of the III and have determined the up-down quark mass difference, Am=7.3 
and 5.2 [MeV] for Pi„ =0.4 and 0.5. 

We have calculated Aa for the CSB parameters fixed by AM. Our results are Aa =0.8 
and 0.5 [fm] for Pni =0.4 and 0.5, which are smaller than the observed value. It is found 

that the contribution of the u-d mass difference to Aa is comparable with that from EM 
interaction because the contributions of OGE and the quark kinetic energy cancel out each 
other. It is pointed out that Aa is sensitive to the choice of the quark model parameters 
because of this cancellation. 
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The P-wave CSB observable, AA{6), is calculated for Pni = 0.4. It is found that CSB of 
the short range part in nuclear force is large enough to explain AA{6). This result is different 
from the conclusion of Brauer et al . We have found that this discrepancy is attributed 
to the introduction of the quark symmetric spin-orbit interaction and the erroneous choice of 
the 7i in their paper. We also have investigated the importance of individual mixing matrix 
element, {^Pi\T\^Pi) , {^D2\T\^D2) and (^F^lTl^Fs) and also the relative importance of the 
OGE and EM interaction. It is found that the contributions of {^Pi\T\'^Pi) and OGE are 
dominant in the observed 6 region. Future experiments for other angles as well as different 
energies may give us further information of the mixings of other partial waves and properties 
the spin-orbit parts of the OGE and EM interactions. In fact, we have observed that at 
En = 477 [MeV] the contributions of the higher partial waves become more important than 
at En = 183 [MeV]. The present quark model description is found to account for the short- 
range part of CSB. We would like to stress that the CSB for the single nucleon as well as 
the central and spin-orbit parts of the nuclear force are consistently described. There is a 
possible remaining short-range contribution introduced by Goldman et al. in Ref.0 (GMS), 
which comes from interference between the QCD and QED effects. GMS pointed out that 
such an interference is necessary to explain the mass difference of the neutoral and charged 
pions. Its effect on the NN scattering was studied by Kao and Yang js^. Because this effect 
has much ambiguity, we have not included its effect in the present sutdy in order to see how 
the current data can be accounted without such complex effects. 

Effects of longer range CSB may require further analysis. Approaches based on the 
chiral effective theory were performed in Refs. |^ . Although the applicability of the chiral 
perturbation theory at high energy NN scattering phenomena is not established, its extension 
to the spin-orbit interaction might be interesting to pursue, which is a subject for future 
works. 

APPENDIX A: THE DECOMPOSITION OF THE T-MATRICES 

The representations of the T-matrices in the basis of the nucleon spins are shown explicitly 
in Appendix A. First we expand the wave function of the two nucleons as 

\p,s:,s',) = Va^Y. ' E'''(L,L„S,S,|J,J,)rs+%)YL,u(p)|.:,4) (Al) 

L,S,J Lz+Sz^Jz 
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Using the wave function Eq. ()A1|) . we calculate the T- matrix. For example, the T-matrix of 
the ^Pq —>-^ Pq scattering is given by 



3Po^='Po 



471 (1, m, 1, s,|0, 0)*(1, 0, 1, 0|0, 0)Yi,^(fc)*Yi,o(p) 
{'Fo\T\'Fo){sl,st\s:,s':'^ - 



-s e"*"^ -s e"*"^ ^ 



c c 
c c 

s e*"^ s e'^ 



V 



(A2) 



where p is the unit vector along the initial momentum p and we take it along the z-axis. We 
show the the T-matrix of each partial in terms of s = sin 6', c = cos 6' and (p, where {6, (p) is 
the scattering angle in the center of mass system. 



Ti 







1 -1 
-1 1 



(A3) 



1 1 
1 1 



(A4) 



1 -1 
-1 1 



(A5) 



3Po^3po 



Tspo 



-s e"'* 


—s e~ 


c 


c 


c 


c 


s e''^ 


s e''^ 



i4> 



(A6) 
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3P2-»3p2 



T3 



3D3^3D3 



Ti 



2c 

s e'"^ -s e-^^ 

s e*'^ -s e^*"^ 

2c 

6c 2s e" 

-3s e^<^ 4c 

-3s e^<^ 4c 



/ 

2s e- 

4c 
4c 







7Id2_^id2 = -TiD2(3c^-l) 



-2s e">' -2s 



1 -1 
-1 1 





3s e-'"t' 
3s e-*"^ 
6c 



V 



^Di-*3Di — tTsDi 



(3c2 - 1) -6sc e-'^ -6sc e''^ ?,s^ e'"'^'*' 



3sce^<^ 2(3c2-l) 2(3c2 - 1) -3sc e"" 
3sce^^ 2(3c2-l) 2(3c2 - 1) -3sc e"^^ 
3s2 e^'-^ 6sc e''^ 6sc e^"^ (Sc^ - 1) 

^ (3c2 - 1) -s2 e-*2,^ ^ 

sc e*'^ -sc e^^'^ 

sc -sc e"^"^ 
-s^e^^-^ (3c2-l) y 

4(3c2 - 1) 6sc e-'"^ 6sc e''^ 2s^ e"'^' 

-8sce^^ 3(3c2-l) 3(3c2 - 1) 8sc e"^ 

-Ssce^"^ 3(3c2-l) 3(3c2-l) 8sc e"^ 

^ 2s^ e^''^ -6sc e''^ -6sc e^"^ 4(3c2 - 



^TiF3(5c3-3c) 







1 -1 
-1 1 
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V 



c(5c2 - 3) 
s(5c2 - 1) e''^ |c(5c2 - 3) 
s(5c2 - 1) e 



-|s(5c2 - 1) e-^"^ -is(5c2 - 1) e' 



lc{5c^ - 3) 



|s(5c2 - 1) 



[c(5c2 - 3) 
\c{5c^ - 3) 



|s(5c2 - 1) e^' 



/ ic(5c^-3) 



;T3F3 



|s(5c2 - 1) e^-^ -|s(5c2 - 1) e"^"^ 

|s(5c2 - 1) e^^ -|s(5c2 - 1) e"*"^ 



lc{5c' - 3) 



-s(5c2 - 1) e-''t' 
-s(5c2 - 1) e"^'^ 
c(5c2 - 3) y 
(A14) 

(A15) 



V2 



3c — 1 3sc e" 



3sc e' 



2 „-i(b \ 



3s' e 



V2^ 

— T3D^_>3Si 



V6 



T3p2^3F2 



3^6 



TlPi_>3pi 



3sce''f' -(3c2 


- 1) -(3c^ - 1) - 


3sce'^ -(3c2 


-1) -(3c2-l) - 


3^2 g2# _3g^ g# _ 


3sc e^'t' 


'i \ 






-1 -1 






-1 -1 












c(5c2 - 3) 


s(5c2 - 


1) e"*^ s( 


s(5c2 - 1) e^"^ 


-c(5c2 


-3) 


s(5c2 - 1) e^"^ 


-c(5c2 


-3) 




-s(5c2 - 


- 1) e'"^ - 


2c — s e"* 


</> _5 g-» 


^ 


-s e^'^ -2c 


-2c 


s e-''t' 


-s e^^ -2d 


-2c 


s e'^'f' 


s e*"^ 


s e*"^ 


2c ^ 


^ -s e-''^ s 






















^ -s e'^ se''^ ^ 





(A16) 



-c(5c2 - 3) 
-c(5c2 - 3) 
s(5c2 - 1) e 



i4> 



(A17) 

-s(5c2 - 1) e-^"^ 
-s(5c2 - 1) e-^"^ 
c(5c2 - 3) I 

(A18) 

(A19) 



(A20) 
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Tz 



3^6 





s e'^ s e-'^ 

-s e*"^ -s e-''^ 




^ -sc e-'^ sc e-'^ ^ 



Ti 



5^6 



T3 



5^6 



T3D2^1D2 



r - 







-sc e*"^ sc e'*" 



sc e^'^ sc e-''^ 

-sc e'"^ -sc e"^' 


-s * (5c2 - 1) e-^<^ s * (5c2 - 1) e'^'*' 


\^ -s* (5c2 - 1) e^'f' s * (5c2 - 1) e'^ j 

^ 

s(5c2 - 1) e^'t' s(5c2 - 1) e"*"^ 

-s(5c2 - 1) e*<^ -s(5c2 - 1) e'"' 




(A21) 



(A22) 



(A23) 



(A24) 



(A25) 



Substituting the above T- matrices into the denominator and the numerator of Eq. (j59p . we 
obtain 

Tr[TtT] = ^|-2TiSo + 2T3s, -2v^T3s,^3D, 

+ 2c(— 3Tip^ + T3pp + 2T3p2 — V6T3p2^3p2) 

+ (3c2 - l)(-5TiD, + 2T3D, + 3T3D3 - V^Tag^^aoJ 

+ c(5c2 - 3)(-7Tip3 + 3T3P, - V6T3p,^3pJ|2 
+ o l^'^^So + 2T3Si - 2V2T3Si_^3Di 

o 

+ 2c(3Tip^ + T3p(, + 2T3p2 — -\/6T3p2_^3p2) 

+ (3c2 - l)(5TiD2 + 2T3D, + 3T3D3 - V2T3s,_3dJ 
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+ c(5c2 - 3)(7TiF3 + 3T3F, - VeTap^^spJI' 

+ i|4T3si + 2V2T3s,^3D, + 2c(3T3p^ + 3T3p, + y6T3p,^3Fj 

o 

+ (3c2 - l)(T3Di + 5T3D2 + 4T3D3 + V2T3s,^3dJ 
+ c(5c2 - 3)(2T3F, + 7T3F3 + v^T3p,_,3fJ|' 
+ ^s^|3T3D^ - 5T3D, + 2T3D3 + 3\/2T3s,^3D, 

o 

35 

+ c(10T3F2 - yTsFj + 5V6T3p,^3fJ|' 

+ ^S^|2T3p„ - 2T3P2 + V6T3P2^3F2 

+ 3c(2T3Di - 2T3D3 - V2T3s,_3dJ 
+ (5c2-l)(3T3F2- v^T3p,^3fJ|' 

+ ^S^|3T3p^ - 3T3P2 - V6T3p2^3p2 

+ c(3T3Di + 5T3D2 - 8T3D3 + 3-\/2T3Si-^3dJ 

+ (5c' - 1)(2T3F2 + ^T3F3 + v^T3p2_3fJ|' (A26) 

1 1 /s 

Tr[Tt((7n - (7p)ATcsb] = -|i(3V6sT3p,_,ip, + bV^scTz^^^i^, + ^s(5c' - 1)T3f3-.if3) 

{4Tis„+4T3s,+2V2T3s,^3D, 

+ 2c(3T3p^ + 3T3P2 + v^T3P2^3fJ 

+ (3c2 - l)(10TiD2 + T3Di + 5T3D2 + 4T3D3 + 3V2T3Si_3di 

+ (5c' - 3)c(14TiF3 + 2T3F2 + 7T3F3 + ^T3P2^3fJ 

+ s'(3T3D^ - 5T3D2 + 2T3D3 + 3a/2T3Si^3dJ 
35 

+ S'C(10T3F2 - yT3F3 + 5V6T3p2^3F2)} (A27) 
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